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Abstract 

We investigate the effect of vector quarks on the inclusive decays B — > X^sl- We show 
that the branching ratio of B — > can differ sizably from the SM and MSSM pre- 
dictions, being enhanced to present experimental observability or suppressed such that 
present runs of the B factories would not observe it. Current measurements of the direct 
CP asymmetry (Acp) for B —* X s ^f are sensitive to the contribution from B — » X^-y. 
For a sufficiently enhanced BR(S — > X^j) we show that the dominant contribution to 
the combined asymmetry may be from B — > X^j. Thus any large value for Acp should 
not immediately be attributed to B — > X s ^, which stresses the importance of good K/ir 
separation. 
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1 Introduction 



Theoretical studies of rare decays of b quarks have attracted increasing attention since the start 
of the physics program at the B factories at KEK and SLAC. Both B factories are running to 
expectations and in excess of 30 fb _1 of data has been accumulated by each experiment. The 
much anticipated measurement of sin20i has established CP violation in the B system 

Many rare decays will be observed for the first time over the next few years, and theoretical 
studies of such decays in the context of models with new physics will continue to intensify. In 
this paper we are concerned with the decays b — > d'y and b — > S7 in a model with vector (or 
singlet) quarks. In such a framework the CKM matrix is necessarily non-unitary, leading to 
flavour changing neutral currents at tree-level. 

It has been known for a considerable time that the inclusive decay B —>■ X s 7 is a sensitive 
probe of new physics @. It has been measured at CLEO [gg, ALEPH [g and BELLE 0. The 
current branching ratio (BR) is in agreement with the Standard Model (SM) expectation M but 
leaves room for new physics contributions. Preliminary measurements of the CP asymmetry 
have been made in the inclusive channel B — > || and the exclusive channel B — > K*^ [[LT]] , 
fLlf . Results with higher precision are expected to come from the B factories in the future, 
especially after a possible luminosity upgrade. 

Transitions of the form b — > d^ have so far remained unobserved and in the SM are sup- 
pressed relative to b — > 57 by a factor |T4d/V^| 2 « 1/20. Experimental upper limits exist for the 
branching ratios (BRs) of the exclusive decay channels, B — > p°7 and B — > p +r ). CLEO |L0] 
obtains < 1.7 x 10~ 5 and < 1.3 x 10 -5 respectively, with corresponding measurements by 
BELLE |T| of < 1.06 x 10~ 5 and < 0.99 x 10" 5 . 

There is considerable motivation for measuring the BR and CP asymmetry (Acp) of the 
inclusive channel B — > X,f)'- 



(i) It provides a theoretically clean way of measuring Vtd, as proposed in |p~2|, |1~3 

(ii) Acp in the SM is sizeable, and much larger than that for b — > 57 ||13| . 



(iii) Acp is sensitive to new physics at the weak scale [|14|-|18 



(iv) b — > d*y transitions sizably affect the measurements of Acp for b — > 57 0. Therefore 
knowledge of Acp for b — >• dj is essential, in order to compare experimental data with 
the theoretical prediction in a given model [|1^ . 



(v) Acp for the combined signal of B — > X s ^ and B — > Xffi is expected to be close to zero 
in the SM [^9,|20|, due the real Wilson coefficients and the unitarity of the CKM matrix. 
Both of these conditions can be relaxed in models beyond the SM. 

The exclusive decays (e.g. B — > p°,u; 7 and B — > p +r y) are expected to be observed for 
the first time at the B factories. A measurement of the inclusive decay, although challenging, 
may be feasible due to the improved K/tc separation which is necessary to reduce the large 
background from b —> sj decays. The much larger Acp of b —>■ d'-f with respect to b — > 57 
(in the SM) is expected to compensate for its smaller branching ratio in terms of experimental 
observability. 



The contribution of vector (i.e. 5C/(2)-singlet) quarks |2"T|, |2"2] to B physics processes has 
received increasing attention over the last few years [p3H30fl. Such particles appear in Grand 



Unified Theories (GUTs) p3[ and mix with the ordinary quarks, thus rendering the CKM 
matrix non-unitary. We will be working in the context of a model with one U- and one D— 
type vector quark. The presence of such particles affects the decays B — > X s j and B — > 
X d j. Due to less stringent experimental constraints on the VQ parameters, the latter decay 
in particular may have a BR and Acp very different from that expected in the SM and the 
Minimal Supersymmetric Standard Model (MSSM). 

Our work is organised as follows: In section 2 we introduce our formalism and the vec- 
tor quark model (VQ-model). In section 3 we outline our approach to calculate the VQ- 
contributions to the BRs and Acp of b — > dj and b — > sj. Section 4 presents the numerical 
results and section 5 contains our conclusions. 



2 Vector quarks and the decays b — ► s,d~f 

There is much theoretical and experimental motivation to study the ratio 

BR(B -> X dl ) 

BR(B - X sl ) {L> 

because it provides a clean handle on the ratio |Vtd/V^| 2 ||13|| . In the context of the SM, R is 
expected to be in the range 0.017 < R < 0.074, corresponding to BR(£> — > X d ^) of order 10~ 5 . 
R stays confined to this range in many popular models beyond the SM. This is because new 
particles such as charginos and charged Higgs bosons in the MSSM contribute to b — > s(d)j 
with the same CKM factors. Therefore C-j is universal to both decays and cancels out in the 
ratio R. In a model with vector quarks this is not the case, and we shall see that R can be 
suppressed or enhanced with respect to the SMQ 

The dominant source of CP-Asymmetry in the decay modes B — > X s /cTj is direct CP- 
violation. Mixing induced CP-violation is strongly suppressed by a factor of m s / d /mi ) . The 
interference between mixing and decay necessary for this kind of CP-violation can only occur 
between identical final states. The photons from b — > s/dj are predominantly left-handed, 
while those from b — ► s/d'y are predominantly right-handed, so interference is possible only 
between strongly suppressed contributions. (Since the weak interaction is left-handed, the spin 
flip in the quark line in the penguin loop must occur on an external leg, and this is proportional 
to the mass of the external quark involved.) 

Mixing induced CP-violation becomes important for B — > X s / d , ~f only if new, right-handed 
weak interactions are possible, because then the spin flip can occur inside the loop EH] . This is 



not relevant for the model we consider because our Vector Quarks are 5'L r (2)-singlets and do 
not introduce right-handed couplings to the weak bosons. 
The direct CP-Asymmetry is given by 

A Ms,) = r(B^X dia) 7)-T(B^X 1 M>y) = AT d(s) 
F(B - X d{s)1 ) + T(B - X-mo) Tf 



1 SUSY models with non-flavour diagonal SUSY-breaking terms can also suppress or enhance R |32,p 



2 



In the SM Af7 p is expected to lie in the range —7% < Ap P < —35% [Tj|, where the uncertainty 



arises from varying the Wolfenstein parameters p and r\ in their allowed ranges. Also included 
is the dependence of AJ P on the scale //& which arises from varying m^jl < [L\> < Im^. For 
definiteness we fix fi b = 4.8 GeV, and find —5% < Af7 p < —28%. Therefore Ap P is much larger 
than A S J P (< 0.6%). 

If b — > c?7 and 6 — > 57 cannot be properly separated, then only of a combined sample 
can be measured.^] Neglecting the masses m s and and sub leading CKM-factors V* s ^V u b 
against V t *^V t b (i.e. considering only the leading terms), it has been shown |T9|,|2"0"f that A S J P 

and Ap p cancel each otheiQ. 

In the presence of new physics such a cancellation does not occur, as was shown in [18 in 
the context of the effective SUSY model. As stressed in JTJ], a reliable prediction of A d Q P in a 



given model is necessary since it contributes to the measurement of A S J P . The CLEO result M 
is sensitive to a weighted sum of CP asymmetries, given by: 

A™* = 0.965^ 7 P + 0.02A d J P (3) 

The latest measurement stands at —27% < A e ^ P < 10% (90% C.L.) ||. The small coefficient 
of Ap P is caused by the smaller BR(5 — > Xjy) (assumed to be 1/20 that of BR(5 — > X a j)) 
and inferior detection efficiencies. 

If the detection efficiencies for both decays were identical, this measured quantity would 
coincide with the weighted sum of the asymmetries 

^s-y+dy _ BR 'jAJp + BR 7 ^4 c 7 p 

cp - BR S7 + BR d7 ' 1 J 

The two terms in eqs. (HH) can be of equal or of opposite sign, i.e. they can contribute 
constructively or destructively to the combined asymmetry. Since BR(I? — > X^) may be 
enhanced in the VQ-model, the relative strength of its contribution to A s ^p dl may be increased. 
The non-negligible contribution of b — > d'y to this combined asymmetry should be verifiable 
at proposed future high luminosity runs of B factories. For integrated luminosities of 200 fb _1 
(2500 fb _1 ), anticipates a precision of 3% (1%) in the measurement of A e ^p. 



We shall be working in a model with an extra generation of vector quarks, U and D. Both 
U and D are singlets under SU (2) and in the interaction basis (denoted by U' and D') do not 
couple to the W . Their mass is generated by terms of the form: 

- fffiM - fuW L U$ + h.c. , (5) 

where f l4 (i = 1— >3) are Yukawa couplings and ip* L = (Ui, d~i)L- Thus U' (D 1 ) will mix with the 
up (down) type quarks, resulting in the (undiscovered) mass eigenstates U and D. The known 
quarks hence contain small amounts of the U' and D' weak eigenstates. A feature of the vector 

2 This is also true (though much less important numerically) for the total rate, where the b — > dry contribution 
should be subtracted from the "6 — > S7" sample as done e.g. in || 

3 An improved analysis with non-vanishing quark masses showed that this cancellation still holds to a very 
high degree for both the inclusive and exclusive decays Bq] 
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quark model is the extended 4x4 CKM matrix which is now not unitary. This can be clearly 
seen from the charged current interaction in the interaction basis: 



C 



V2 



(6) 



whereQ 



J 7 ' = u' L a w ^d' L = wlVckm7^l and u L 



( u L \ 
\UlJ 



( d L \ 

sl 
b L 
\D L 



(7) 



Here the matrix aw is Diag(l, 1,1,0), reflecting the fact that the charged current couplings 
to the known quarks are flavour diagonal, while the vertex W-U'-D' is absent. In the mass 
basis the vertices W-U-d (d denoting any down-type quark) are generated and thus an extra 
diagram with an internal U quark contributes to the decay b — > sj, dj. The 4x4 CKM matrix 
is then simply given by 

(8) 



where U 



ut 



Vckm — U^awUt 
and Ul rotate the mass eigenstates to the interaction eigenstates. 



Clearly V C km 

is not unitary due to the presence of aw- A consequence of this is the generation of FCNC 
vertices, where the effective coupling for the vertex Z-di-dj is (V*V)ij. This gives rise to a 
further class of diagrams that contributes to b — > sj, dj (c.f. section 3). 

These FCNC vertices are strongly constrained by current experiments. The vertices Z-b- 
s and Z-b-d are e.g. constrained by the non-observance of B — > X s l + l~ and B — > Xdl + l~, 
respectively ||25|| . (Note that the first evidence for B — > Kl^^T has recently been reported by 
the BELLE collaboration j|55|.) 

We now briefly summarise the previous works which considered the effect of vector quarks 
on B —>■ X(i(s)l- (Work on these decays in other models has already been summarised in ||18|| .) 



2q,|27l studied a model with a D type vector quark which induces Z and h FCNC con- 
tributions to B — > Xd( s )^- It was shown that these contributions have negligible impact on 
B — > Xg'-f, but can enhance or suppress R. Implications for A S J P were found to be small. 

The effect of the U-W contribution on A S J P was considered in [EC]. In a model with both 
U and D vector quarks (where the D-contribution could actually be neglected) [p8|| , the U- 
W mediated contribution to BR,(B — > X^) was shown to permit BRs anywhere within the 
experimental limits. More accurate measurements of BR(5 — > X s 7) at the B factories will 
further restrict the available parameter space for VQ-models. 

To our knowledge the combined effect of U and D vector quarks on BR(I? — > X<p() has not 
yet been considered. In addition, no analysis of A d ^ P in any model with vector quarks has been 
carried out. Given the importance of the measurement of the combined asymmetry A s ( ]p d ' y , we 
will also study the correlation of the individual asymmetries A S J P and A d ^ P as a function of the 
VQ parameters. 



It is obvious where ul denotes the 4-vector of up-type quarks and where it denotes the up-quark only 
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Z,H 



7 



u, c, t, U 



d, s, b, D 



s, d 



s,d b 
Figure 1: Diagrams contributing to b —>■ s,d^ in the VQ-model 

3 Direct CP Asymmetry in B — > ^d,s1 

In this section we explore the effect of U and .D-type vector quarks on both BR(S — > X dtS l) 
and the direct CP asymmetries, A(fp '. Diagrams contributing to b —>■ s,dj in the VQ-model 
are shown in fig. [[]. 

The effective Hamiltonian for b — > d^ is given by 



K 



* 1=1 



(9) 



where Qi(fib) are the current density operators for the AB = 1 transitions and Ci(fib) are their 
Wilson coefficients. The relevant operators for b —>■ dj decay are given by 



Q2 = d L ^c L c L ^b L , 

g _ 
Q7 = Y^rnbdLa^baF^, 

Qs = ^m b d L a^T a b R G%. 



(10) 



The analogous formulae for the b — > sj decay can obtained from Eqs.(|9|) and ( JTOj ) by making 
the replacement d — > s. 

We use formula (27) from pO|j : 



82 

xlm[(l + e dW + A d{s) )C 2 C* 7 ] + — b(z, 5) Im[(l + e d{s) + A d(s) )C 2 C* 8 



(11) 



where in our notation = (z x b — V^ x Vub) / (V* x Vtb) for (x = d,s). 

The expressions for the Wilson coefficients at the Mw scale may be divided into the charged 
current and neutral current mediated contributions, where only the former exist in the SM. 
Contrary to the SM case, the charged current mediated contribution contains a term stemming 
from a loop with the U vector quark. 

Thus at the Mw scale we have: 



C 7 = + C$ 



iZH 



(12) 
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We will use the leading order expressions for CV (with the constant terms included in the 
Inami-Lim functions) which can be found in |29| . 



In the SM the CKM-factors are customarily included in the pre-factor of the effective 
Hamiltonian. This is possible because the top-quark loop is the dominant contribution to C7 
and the CKM-factor V* x V c b in C 2 can be reexpressed in terms of V t * x V t b utilizing CKM unitarity. 

In the vector quark model this is not possible anymore. For CV, the constant terms in the 
Inami-Lim functions that are cancelled in the SM due to CKM unitarity, have to be taken into 
account. Also in replacing V* x V c b in C2, additional terms arising from CKM non-unitarity are 
generated. This results in different Wilson coefficients for b — > 37 and b — > d'y. 

For comparison with standard calculations, we keep the term V t * d ^Vtb in the pre-factor of 
Heff and therefore have to divide all but the top-quark contributions to the Wilson coefficients 
by this factor. 

We define Cp and Cp as the Wilson coefficient for b — > 57 and b — > d'y respectively. The 
magnitude of \Cp\ is constrained by measurements of BK(B — > Xg'y) ||28|| . 

In contrast, Cy 7 is only very weakly constrained due to the non-observation of b — > d'y. 
Recent direct searches for the exclusive decay B — > fry give an improved upper bound on -R ex ci 



for exclusive decays of -R ex ci < 0.19 (90% C.L.) [ill]. One expects a weaker bound on the 
inclusive R, since estimates of R exc \/R are smaller than 1 ||27|| . 

On the other hand the presence of b — > dj events in the samples of the inclusive measure- 
ments of BR(5 — > Xg'y) || also constrains R. This bound, however, is not relevant for our 
analysis, since in the vector quark model the prediction for b — > 57 can easily be lowered to 
account for a higher admixture of b — > d'y in the experimental sample. 



4 Numerical Results 



We vary the vector quark parameters in the ranges given in table |l| The Wolfenstein parameters 
p and r\ quantify the current uncertainty in the CKM parameters related to the mixing of the 
first and the third generation of quarks. Our results are virtually independent of our choice of 
rriu/D and M H . 

In the s-sector we use |V^ a Vt/&| < 0.004 which is 1/3 smaller than the bound derived 
from BR S7 in [28], we take this smaller value in the light of recently improved measurements 
The non-unitarity parameter z s b {zdb) is constrained from non-observation of 
(B^X d £+£-) [||. 

Our results depend crucially on the magnitude of |V^ d Vj/b|. Bounds from CKM unitarity on 
this quantity are rather weak (of the order of the uncertainty in the CKM parameters in the d- 
sector) and B — ► X^j has not yet been observed, so the only restrictions on this parameter come 
from i?°i? -mixine' The exact bounds on \ V T *,Vnh\ from B°B° -mixing are not clear and 



of b —>■ 57. 
B -> X s £ + £ 



scing p2j-p4|. The exact bounds on | V^ rf Vf/b| 
depend on the tree level Z-mediated contributions as well as the rather large uncertainties in 
the SM prediction of B°B° -mixing (bag parameter, etc.). We give results for different choices 
of |V^ d V[/&| up to 0.01 (« |V^Vtf,|). A more detailed study of 5°_B°-mixing in the VQ model 
will be presented elsewhere J57 . 



In figure @ we plot Ap P against Aq P . It can be seen that while A h J P does not substantially 



differ from its SM value, A^p can vary over a much larger range. We restrict ourselves to 
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-0.25% 0.25% 0.5% 0.75% 1% 1.25% 

Figure 2: CP asymmetry of b — > d'y against CP asymmetry of b — > sj. The VQ-parameters 
are varied in the ranges given in table |]; each point in the diagram corresponds to one random 
point in parameter space. 
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Table 1: Model parameter ranges 
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Figure 3: Contour plot of CP asymmetry against branching rate for b — > d'y. Contours for 
different values of |V^ d V[/b| are given, and Arg V^ d Vub is varied along the contours 



points in parameter space where I-A^pI < 45%. Asymmetries greater than 50% are attainable, 
but these correspond to cases with a virtually complete cancellation of the SM contribution by 
the new VQ contribution (IC^ 7 ! <C 0.3). Such large asymmetries are untrustworthy, since our 
formula for the CP asymmetry (Q) breaks down in these cases. 

The correlation between Aq P and BR d7 is studied in detail in figure |3|, where it can be seen 
that \A%\ > 45% occurs only for BR d7 < 10 6 . Branching ratios of this magnitude would 
require 3> 10 8 bb pairs to be detected which is beyond the discovery potential of current B 
factories. Note that figure [3] was obtained with fixed values for the CKM parameters p and rj 
as given in table |I| Varying p and rj shifts the contours in figure [5] vertically by the amount of 
uncertainty in the SM prediction of A^'p. 

In figure |] we plot the ratio of the branching ratios for b — > d'y and b —>■ sj against 
the argument of V^ d Vub, drawing curves for different values of |V^ d Vj/&|. For real, positive 
Vy d Vubi there is constructive interference with the SM contribution proportional to V t * d V t b = 
AX 3 (1 — p — irj) + 0{\ 4 ). For real, negative V^ d Vub (corresponding to points around Arg 
^YudYub) = 7r) the interference is destructive. Note that the smallest values for R are not 
obtained for the largest |V^ d Vt/&|, because if |V^ d V[/f,| > |V^V 4 b|, exact cancellation cannot 
occur. 

In figure || we plot the combined CP asymmetry as defined in equation (H) against the 
argument of V^ d Vub- We plot several curves for different values of |V^ d V[/&| just like in figures 
and f|. The extreme values for the combined asymmetry are obtained where R is maximal 
and not — as could be naively expected — where the CP asymmetry of b — > d'y becomes much 
larger than its SM value. This is because large CP asymmetries in b — > d'y as seen in figure [3] 
always imply small values for R which makes these points in parameter space unimportant for 
the combined asymmetry. 

Note that in our analysis BR/ 7 and A b ^p sl are close to their SM values. The huge variations 
in A s Jp d,y stem from the variation in BR d7 . In wide ranges of our parameter space, b — » dy 
actually dominates the combined asymmetry ! Any large signal observed in ^4^7p d7 (which is 
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Figure 4: R against Arg V^ d Vub for different values of |V^ d V[/b| 
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Figure 5: Combined Asymmetry against Arg V^ d Vub for different values of (V^V^I 



the experimentally relevant quantity) should not necessarily be attributed to b — > 57 only This 
motivates excellent K-tt separation which is crucial in distinguishing the two decay modes. 

To obtain a rough estimate on the number of bb pairs N bb required to establish CP violation 
to a given significance, following |1| we define the observability f^( s ) 7 = (A 1 ^^ 1 ) 2 ■ BR^-* 7 , 
where fld( s )-y oc 1/N bb . This observability is plotted in figure [| against Arg V^ d Vub, where the 
horizontal lines indicate the SM observabilities. (We do not take into account different detection 
efficiencies for the two channels.) Despite the suppressed BR d7 in the SM, A b ^p dl is much more 
observable than A b ^p sl . In most of the parameter space, A b ^p dl has a higher observability than 
in the SM. The optimal observability occurs for the largest values of R (compare fig. f|). 
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Figure 6: CP asymmetry and observability Q against Arg V^ d Vub 

5 Conclusions 

We have studied the effect of a single generation of vector quarks (VQ) on the short-distance 
component of the inclusive decays B — > Xd tS j. We found that A h ^p sl is relatively insensitive to 
VQ interactions if the experimental bounds on BR S7 are respected. For the decay mode h — > 
both the decay rate and the CP asymmetry can be significantly different from their SM values. 
Phenomenologically most salient is the fact that the decay rate can be increased up to current 
experimental observability or decreased beyond the sensitivity of the current B factories. This 
variation of BR* 7 has a profound impact on the experimentally relevant CP asymmetry of a 
combined sample of b — > s,dj which can actually be dominated by A 1 ^ 1 . Observation of 
any A s Jp dl significantly different from zero is a clear sign of new physics, but the theoretical 
interpretation should await experimental separation of the individual channels. 
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